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Abstract. This paper deals with the partition function of the Ising model 
from statistical mechanics, which is used to study phase transitions in phys- 
ical systems. A special case of interest is that of the Ising model with con- 
stant energies and external field. One may consider such an Ising system as 
_ a simple graph together with vertex and edge weights. When these weights 

are considered indeterminates, the partition function for the constant case is a 
• trivariate polynomial Z(G;x,y, z). This polynomial was studied with respect 

' to its approximability by L. A. Goldberg, M. Jerrum and M. Paterson in I14| . 

Z(G; x, y, z) generalizes a bivariate polynomial Z{G; t, y), which was studied in 
by D. Andren and K. Markstrom in yj. 
^ ■ We consider the complexity of Z{G; t, y) and Z{G; x, y, z) in comparison to 

that of the Tutte polynomial, which is well-known to be closely related to the 
' Potts model in the absence of an external field. We show that Z(G;x,y,z) is 

#P-hard to evaluate at all points in Q^, except those in an exceptional set of 
low dimension, even when restricted to simple graphs which are bipartite and 
planar. A counting version of the Exponential Time Hypothesis, #ETH, was 
' introduced by H. Dell, T. Husfcldt and M. Wahlen in 6 in order to study the 

, complexity of the Tutte polynomial. In analogy to their results, we give under 

• ■ #ETH a dichotomy theorem stating that evaluations of Z(G;t,y) either take 

^ ' exponential time in the number of vertices of G to compute, or can be done 

in polynomial time. Finally, we give an algorithm for computing Z{G; x, y, z) 
in polynomial time on graphs of bounded clique- width, which is not known in 
^•f^ ■ the case of the Tutte polynomial. 

>: 

m 
rn 

■ An Ising system is a simple graph G — (V, E) together with vertex and edge 
weights. Every edge {u, v) G E has an interaction energy and every vertex u V 
has an external magnetic field strength associated with it. A function a : V {il} 
is a configuration of the system or a spin assignment. The partition function of an 
Ising system is a generating function related to the probability that the system is 

_ in a certain configuration. 

■ L. A. Goldberg, M. Jerrum and M. Paterson 14 studied the Ising polynomial 
d I in three variables Z{G; x, y, z) for the case where both the interaction energies 

of an edge {u, v) and the external magnetic field strength of a vertex v are con- 
stant. They consider the existence of fully polynomial randomized approximation 
schemes (FPRAS) for the graph parameters Z{G;j,S,e), depending on the values 
of (7, 6, e) G Q^. They provide approximation schemes for some regions of Q3 while 
showing that other regions do not admit such approximation schemes. Approxim- 
ation schemes for Z(G;x, y,z) were further studied in [351 [2S]. M. Jerrum and A. 
Sinclair studied in [T^ the approximability and #P-hardness of another case of the 
Ising model, where weights are provided as part of the input and no external field 
is present. The bivariate Ising polynomial Z(G;t,y), which was studied in T for 
its combinatorial properties, is equivalent to setting x = z = t in Z(G;x, y,z). It 
is shown in [I] that Z{G; t, y) encodes the matching polynomial, and is equivalent 



1. Introduction 



This work was partially supported by the Fein foundation and the graduate school of the 
Technion. 



2 



TOMER, KOTEK 



to a bivariate generalization of a graph polynomial introduced by B. L. van der 
Waerden in ^5] . 

The trivariate and bivariate Ising polynomials fall under the general framework 
of partition functions, the complexity of which has been studied extensively start- 
ing with [7] and followed by [1 [H HZl H] . From gZl Theorem 6.1] and imphcitly 
from [13] we get that the complexity of evaluations of the Ising polynomials sat- 
isfies a dichotomy theorem, saying that the graph parameter Z{G\^,5) is either 
polynomial-time computable or ^^P-hard. However, 5 must be positive here. 

The g-state Potts model deals with a similar scenario to the Ising model, except 
that the spins are not restricted to ±1 but instead receive one of q possible values. 
The complexity of the q-state Potts model has attracted considerable attention in 
the literature. The partition function of the Potts model in the case where no 
magnetic field is present is closely related to the Tutte polynomial T{G;x,y). It 
is well-known that for every 7,(5 £ Q, except for points (7,(5) in a finite union of 
algebraic exceptional sets of dimension at most 1, computing the graph parameter 
r(G;7,(5) is #P-hard on multigraphs, see [5]. This holds even when restricted to 
bipartite planar graphs, see [30] and [H]. In contrast, the restriction of the Tutte 
polynomial to the so-called Ising hyperbola, which corresponds to the case of Ising 
model with no external field, is tractable on planar graphs, see [9l ITSl [8]. 

H. Dell, T. Husfeldt and M. Wahlen introduced in [6] a counting version of 
the Exponential Time Hypothesis (#ETH), which roughtly states that counting 
the number of satisfying assignments to a 3CNF formula requires exponential time. 
This hypothesis is implied by the Exponential Time Hypothesis (ETH) for decision 
problems introduced by R. Impagliazzo and R. Paturi in [TO. Under ^j^ETH, the 
authors of [6] show that the computation of the Tutte polynomial on simple graphs 
requires exponential time in j^^^^ in general, where nia in the number of edges of 
the graph. For multigraphs they show that the computation of the Tutte polynomial 
generally requires exponential time in ma- 
in this paper we prove that the bivariate and trivariate Ising polynomials satisfy 
analogs of some complexity results for the Tutte polynomial. For the bivariate Ising 
polynomial we show a dichotomy theorem stating that evaluations of Z{G; t, y) are 
either #P-hard or polynomial-time computable. Moreover, assuming the counting 
version of the Exponential Time Hypothesis, the bivariate Ising polynomial require 
exponential time to compute. Let uq be the number of vertices of G. 

Theorem 1 (Dichotomy theorem for the bivariate Ising polynomial). 

For all (7,5) G Q^' 

(i) If J £ {—1,0,1} or 5 = 0, then Z{G;j,d) is polynomial time computable. 

(ii) Otherwise: 

• Z{G; 7, (5) is ^P-hard on simple graphs, and 

• unless ^ETH fails, requires exponential time in ^^^'^^ on simple 
graphs. 

We show that the evaluations of Z(G;x, y,z), except for those in a small excep- 
tional set B C Q^, are hard to compute even when restricted to simple graphs 
which are both bipartite and planar. 

Theorem 2 (Hardness of the trivariate Ising polynomial). 

There is a set B C such that for every {■j,S,e) € Q^\i3, Z{G;j,S,e) is ^P-hard 

on simple bipartite planar graphs. 

B is a finite union of algebraic sets of dimension 2. 

Although Z(G;x, y,z) is hard to compute in general, its computation on restric- 
ted classes of graphs can be tractable. Computing Z{G; x, y, z) is fixed parameter 
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tractable with respect to tree- width using the general logical framework of ^9] . This 
implies in particular that Z(G;x,y,z) is polynomial time computable on graphs of 
tree- width at most k, for any fixed k, which also follows from [22]. Likewise, the 
Tutte polynomial is known to be polynomial time computable on graphs of bounded 
tree- width, cf. [HHT]. In contrast, for graphs of bounded clique- width, a width 
notion which generalizes tree-width, the best algorithm known for the Tutte poly- 
nomial is subexponential, cf. [12] • We show the following; 

Theorem 3 (Tractability on graphs of bounded clique- width) . 
There exists a function f{k) such that Z{G; x, y, z) is computable on graphs of clique- 
width at most k in running time O (nc^'-'^-') . 

In particular, Z(G';x, y,z) can be computed in polynomial time on graphs of 
cliquc-widtlQ at most fc, for any fixed k. On the other hand, it follows from [TT] 
that, unless FPT = W[l], Z(G;x, y,z) is not fixed parameter tractable with re- 
spect to clique-width^ i.e. there is no algorithm for Z(G;x, y,z) which runs in time 
O {q{nc) • f{k)) on graphs G of clique-width at most k for every k such that g is a 
polynomial. 

2. Preliminaries 

2.1. Definitions of the Ising polynomials. Let G be a simple graph with vertex 
set V{G) and edge set E{G). We denote no = \V{G)\ and ma = \E{G)\. AU graphs 
in this paper are simple and undirected unless otherwise stated. 

Given S C V{G), we denote by Ec{S) the set of edges in the graph induced by 
5" in G and by Eg{S) the set of edges in the graph obtained from G by deleting 
the vertices of S and their incident edges. We may omit the subscript and write 
e.g. E{S) when the graph G is clear from the context. 

Definition 4 (The trivariate Ising polynomial). The trivariate Ising polynomial is 

Z(G;X,y,z)= ^ y,\Ea(S)\^\S\^\EaiS)\ ^ 
SCV(G) 

For every G, Z(G;x, y,z) is a polynomial in Z[x, y,z] with positive coefficients. 

Definition 5 (The bivariate Ising polynomial). The bivariate Ising polynomial is 
obtained from Z{G] x, y, z) by setting x = z = t. In other words, 

Z(G;t,y)- J2 t\^oiS)\+\Eaisys\ _ 

SCV{G) 

The cut [S, S]g is the set of edges with one end-point in S and the other in 
S = V{G) \ S. The bivariate Ising polynomial can be rewritten as follows, using 
that Eg{S), Eg{S) and [S, S]g form a partition of E{G): 

(1) Z(G;t,y)=t™'= t-|[^'^l«lyl^l . 

scy(G) 

The bivariate Ising polynomial is defined in this paper in a way which is slightly 
different from, and yet equivalent to, the way it was defined in [1]. The definition 
in [1] is reminiscent of Equation ([T]) . 

In Section 13] we use a generalization of the bivariate Ising polynomial: 

Definition 6. For every S, G C V{G) such that B n G = we define 

Z(G; B, G; t, y) - ^ ^\Ea(s)\+\Eois)\^\s\ 

BCSCV{G)\C 



"'^ Rank-width can replace clique- width here and in Theorem[3l since the clique- width of a graph 
is bounded by a function of the rank- width of the graph. 
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Clearly, Z(G;0,0;t,y) — Z(G';t,y). In Section H] we use a multivariate version 
of Z{G]B,C;y.,y,2). In Section [S] we use a different multivariate generalization of 
Z(G;x,y,z). 

We denote by [i] the set {1, . . . , «} for every i E 

2.2. Complexity of the Ising polynomial. Here we collect complexity results 
from the literature in order to discuss the complexity of computing, for every graph 
G, the trivariate (bivariate) polynomial Z(G;x, y,z) (Z(G;t,y)). By computing the 
polynomial we mean computing the list of coefficients of monomials x'y-' z*'' such that 
i, G {0, 1, . . . , ma} and j € {0, 1, ... , no}- 

In [T] it is shown that several graph invariants are encoded in Z(G;t,y). 

Proposition 7. The following are polynomial time computable in the presence of 
an oracle to the bivariate polynomial Z(G;t,y). The oracle receives a graph G as 
input and returns the matrix of coefficients of terms t*y^ in Z(G;t,y). 

• the matching polynomial and the number of perfect matchings, 

• the number of maximum cuts, 

and, for regular graphs, 

• the independent set polynomial and the vertex cover polynomial. 

The following propositions apply two hardness results from the literature to 
Z(G;t,y) using Proposition [T] 

Proposition 8. Z(G;t,y) is ffP-hard to compute, even when restricted to simple 
3-regular bipartite planar graphs. 

Proof. The proposition follows from a result in [3T] which states that it is #P-hard 
to compute #3RBP — VC, the number of vertex covers on input graphs restricted 
to be 3-regular, bipartite and planar. □ 

For the next proposition we need the following definition which is introduced in 
[B] following [TB] : 

Definition 9 Exponential Time Hypothesis (#ETH)). Let s be the infimum 
of the set 

{c : there exists an algorithm for =i^3SAT which runs in time 0(c"'^)} 

The ^ Exponential Time Hypothesis is the conjecture that s > 1. 

Proposition 10. There exists c > 1 such that the computation of Z(G; t, y) requires 
n (c"*^ ) time on simple graphs, unless ^ETH fails. 

Proof. The claim follows from a result of [5] which states that there exists c > 1 
for which computing the number of maximum cuts in simple graphs G takes at 
least (c™'^) time, unless the ^^ETH fails. It is easy to see that the problem of 
computing the number of maximum cuts of disconnected graphs can be reduced to 
that of connected graphs and so no subexponential algorithm exists for connected 
graphs, and the proposition follows since for connected graphs no = 0{mG)- 

□ 

On the other hand, Z{G] t, y) and Z{G] x, y, z) can be computed naively in time 
which is exponential in uq. 

The three above propositions apply to Z(G;x, y,z) as well. 
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Figure 1. The graph and the construction of L5 from ^4. 

is obtained from £4 by adding the vertex hd(5) and its incident 
edges, and renaming tri(4) and tr2(4) to tri(5) and tr2(5) respect- 
ively. 



tr2(5) to (4) 




tri(5) t;ri(4) 
(a) (b) 



2.3. Clique-width. Let [k] {1, . . . , fc}. A fc-graph is a tuple (G, c) which consists 
of a simple graph G together with labels G [k] for every v E V{G). The class 
CW{k) of fc-graphs of clique-width at most k is defined inductively. Singletons 
belong to CW{k), and CW{k) is closed under disjoint union U and two other 
operations, Pi^j and iiij, to be defined next. For any i,j e [fc], pi^j{G,c) is 
obtained by relabeling any vertex with label i to label j. For any i, j € [k], fj,i,j{G, c) 
is obtained by adding all possible edges {u,v) such that Cu — i and Ci, = j. The 
clique-width of a graph G is the minimal k such that there exists a labeling c for 
which (G, c) belongs to CW{k). We denote the clique-width of G by cw{G). 

A fc-expression is a term t which consists of singletons, disjoint unions U, rela- 
beling pi^j and edge creations pij, which witnesses that the graph val{t) obtained 
by performing the operations on the singletons is of clique-width at most k. Every 
graph of tree-width at most k is of clique-width at most 2'^+^ + 1, cf. [5]. While 
computing the clique- width of a graph is NP-hard, S. Oum and P. Seymour showed 
that given a graph of clique-width k, finding a (2'^*''+^ — l)-expression is fixed para- 
meter tractable with clique-width as parameter, cf. p51 124) . 

3. Exponential Time Lower Bound 

In this section we prove that in general the evaluations (7, 5) G of Z(G;t, y) 
require exponential time to compute under #ETH. In analogy with the use of 
Theta graphs to deal with the complexity of the Tutte polynomial, we define Phi 
graphs and use them to interpolate the indeterminate t in Z{G; t, y). We interpolate 
y by a simple construction. 

3.1. Phi graphs. Our goal in this subsection is to define Phi graphs and 
compute the bivariate Ising polynomial at y = 1 on graphs G ® to be defined 
below. In order to define Phi graphs we must first define i^-graphs. For every 
h G N, the graph Lh is obtained from the path Ph+i with h edges as follows. Let 
hd{h) denote one of the end-points of Ph+i- Let tii{h) and tr2(/i) be two new 
vertices. is obtained from Ph+i by adding edges to make both tri(ft,) and tr2(/i) 
adjacent to all the vertices of Ph+i- 

We can also construct recursively from L^-i by 

• adding a new vertex hd(/i) to L^-i, 

• renaming tri(/i — 1) to tTi{h) for i = 1,2, and 

• adding three edges to make hd{h) adjacent to hd(/i— 1), tri(/i) and tr2{h). 
Figure [1] shows L5. 

Let U G be a partition of the set {tri, tr2, hd}. Let B(h) be the subset of 
{tri(/i), tr2(/i), hd(/i)} which corresponds to B and let G{h) be defined similarly. 
We have that B(h) and C{h) form a partition of {tri(/i), tr2(/i), hd(ft,)}. 
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Definition 11. We denote bB,c{h) = Z{Lh:, B{h),C{h);t, 1). 
The next two lemmas are devoted to computing hB.c{h). 



Lemma 12. 



^{tri,hd},{tr2}(^) 
^{tri},{tr2,hd}(^) 



&{tr2,hd},{tri}(^) — 
6{tr2},{tri.hd}W = (t^+t)'' -t. 



Proof. We have 



^{tri,hd},{tr2}(^) 
^{tri},{tr2,hd}(^) 



&{tr2,hd},{tri}(^) 
^{tr2}.{tri,hd}(^) 



by symmetry. We compute ^{tri,hd},{tr9}('i) by finding a simple linear recurrence 
relation which it satisfies and solving it. We divide the sum fe{tri,hd}.{tr2}(''') i'^to 
two sums, 

0{tri,hd},{tr2} 

{h) = Z(L;,;{tri(/i),hd(/i),hd(/i-l)},{tr2(/i)};t,l) + 
ZiLh; {tri(/i), hd(/i)}, {tr2(/i), hd{h - 1)}; t, 1) 

depending on whether hd(/i — 1) is in the iteration variable S of the sum 
^{tii,hd}.{tr2}(^) Definition These two sums can be obtained from 

&{tri,hd},{tr2}(^ - 1) and 6{tri},{tr2,hd}(^ " 1) by adjusting for the addition of hd{h) 
and its incident edges: 

• hd(ft, — 1) e 5*: adding hd(/i) (to the graph and to S) puts two new edges 
in E{S) U E{S), namely (tri,hd(/i)) and {hd{h - l),hd(/i)). Hence, 

ZiU; {tnih),Mih)M{h - 1)}, {tr2(/i)}; t, 1) 



• hd(/i ~ 1) ^ S: adding hd(/i) puts just one new edge in E{S) U E{S), 
namely (tri, hd(/i)). Hence, 

Z{Lh; {tri{h)Mih)}, {tr2ih)Mih - 1)}; t, 1) 

= &{tri},{tr2did}(/i - 1) -t- 

Using that 6{tri,hd},{tr2}(^ - 1) = ^{tri},{tr2.hd}(/» " 1), wc get: 

(2) 6{tri,hd},{tr2}(^) = fe{tri,hd},{tr2}(/i " 1) ' (t^ + t) 

and the lemma follows since 6{tri,hd}.{tr2} (0) — t (note Lq is simply a path of length 



We are left with two distinct cases of bB,c{h) to compute, since by symmetry. 



&{tri,hd},{tr2}(^ - 1) • t^ . 



3). 



□ 



O{tri,tr2,hd},0 

{h) 



^0,{tri,tr2,hd}(^) and 6{tri,tr2},{hd}(^) = &{hd},{tri ,tr2 } (^) • 



Lemma 13. Let 




Cl = t^ - C2 , 



t (-t''' - 2 + t + tV5 - 2t2 + t4) 



2V5 - 2t2 + t4 

di = 1 - C?2 , 



-1 - 2t + t2 + V5 - 2t2 + t4 



2V5 - 2t2 + t4 
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Figure 2. An example of a Phi graph: the graph for "H = {1, 3, 4}. 




triCH) 



Ai corresponds to the + case, //t G M then ci, C2, di, (i2, Ai, A2 G M, Ai =/= X2, and 

&{tri,tr2,hd},0(^) = CiX'l + C2X2 
^{tri,tr2},{hd}(/») = diX'l + d2X2 ■ 

Proof. The content of the square root is always strictly positive for t G M. Hence, 
Ai 7^ A2 and ci, C2, c?i, ^2, Ai, A2 G M. 

The sequences 6{tri,tr2,hd}.0(^) and 6{tri,tr2},{hd} (^) satisfy a mutual linear re- 
currence as follows: 

^{tri,tr2,hd}j(/j) = ^{tri,tr2,hd}j(^ - 1) • t^ + 6{tri,tr2}.{hd} - 1) • t^ 
^{tri,tr2},{hd}(/i) = fe{tri,tr2,hd}j(^ ^ 1) + ^{tn ,tr2 } ,{hd} (/i - 1) ' t 

This implies that both fo{tri,tr2,hd},0(^) and fe{tri.tr2}.{hd}(^) satisfy linear recurrence 
relations with the following initial conditions: 

^{tri,tr2,hd},0(O) = t^ and 6{tri,tr2,hd},0(l) ^ + t^ 
fe{tri,tr2},{hd}(0) = 1 and fc{tri,tr2},{hd}(l) ^ t^ + t . 

These recurrences can be calculated and solved using standard methods, see e.g. 
[lO] or [H]. □ 

Using the previous two lemmas, we get: 

Lemma 14. 

Z(L,,;{tri},{tr2};t,l) - {tr2}, {tri}; t, 1) = 2t(t2 + t)'' 

Z{Lh; {tri, tr2}, 0; t, 1) = Z{Lh; 0, {tri, tr2}; t, 1) = (ci + di)A^ + (c2 + d2)A^ , 
where Ci, C2, c?i, ^2, Ai, A2 are as in Lemma \TSl 

Proof. The lemma follows from Lemmas [T^] and [T^l □ 

Definition 15 (Phi graphs). Let H he a, finite set of positive integers. We denote 
by the graph obtained from the disjoint union of the graphs : ft, G "H as 
follows. For each i — 1,2, the vertices tri(/i), h G H, are identified into one vertex 
denoted tr.i('H). 

The number of vertices in <i>^ is 2 + '^Zheu^^ + Figure [2] shows ${1,3,4}. 
Lemma 16. LefH be a finite set of positive integers. Then 

Z{^n\ {tri(H)}, {tr2(H)}; t, 1) = (2t)l«l [] (t^ + t)'' , 

hen 
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and 

Z($„;{tri(H),tr2(H)},0;t,l) = 

^($„;0,{tri(H),tr2(H)};t,l) = X{ {{ci + di)\\ + {c2 + d2)\'i) 



hen 



Proof. It follows from Lemma [fusing that all edges are contained in some Lh- □ 
We can now define the graphs G (^H: 

Definition 17 {G^Ti). Let H be a finite set of positive integers. Let G be a graph. 
For every edge e — (ui, U2) G E{G), let $-^,6 be a new copy of where we denote 
tri(?^) and tr2('H) for ^H,e by tri(-H,e) and tr2('H,e). Let G «) = G «) be 
the graph obtained from the disjoint union of the graphs 

: e e E{G) 

by identifying tr,;('H, e) with m^, i = 1, 2, for every edge e = (^1,^2) G £;(G). 

Lemma 18. Let % he a finite set of positive integers. Let ft^-n and gp^u be the 
following functions: 

/t,-H(ei,e2,ri,r2) = J| (eir^ + e2r2 ) 
hen 



\ hen ) 



Then 

Proof. Let G = G (g) By definition, 

Z(G;t,l)= tl^o(^)^^<5(^)l . 

scy(G) 

We can rewrite this sum as 

Z(G;t,l) = E I n ^(3'w.e;{tri(H,e)},{tr2(H,e)};t,l) 

SCy(G) \ee[5,S]G 

n Z{^n,e\ {tri(H, e), tr2(H, e)}, 0; t, 1) 

\eeEG{S)UEG(S) 

since edges only occur within some ^u^- Using lemma [TOl the sum in the last 
equation can be written as 

\\SAg\ 



Y f(2t)i«in(t^+t)'^) 

scv{G) \ hen / 



n ((ci+di)A'i' + (c2 + d2)A 



\EG{S)uEGiS)\ 



\hen 



^It does not matter how we identify ui and U2 with tri{'H, e) and tV2{'H,e) since the two 
possible alignments will give raise to isomorphic graphs. 
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Since = niG — \Eg{S) U Eg{S)\, we can rewrite the last equation as 



TOG 



(2t)i«i n(t'+t)'j • 

l,,g„((ci+di)A^ + (c2+d2)A5) 




(2t)l«ln,eH(t'+t)' 
The last sum can be rewritten as 



E 

scy(G) 



2t W+tJ 2t Vt^+t 



\EG{S)UEGiS)\ 



\EG{S)UEGiS)\ 



then 

and the lemma follows. □ 

The construction described above will be useful to deal with the evaluation of 
Z(G;t, y) with y = — 1 due to the following lemma. For a graph G, let be the 
graph obtained from G by adding, for each v (E V{G), a new vertex v' and an edge 
{v,v'). So v' is adjacent to v only. G(i) is a graph with 2nG vertices. 

Lemma 19. Z(G;t, 1) = (t - l)-"«Z(G(i); t, -1) 

Proof. By definition we have 

Z(G(i);t,-l) = tl^'=a)(^)^^'=a)(^)l (-1)1^-1 

Scy(G(i)) 

^ i\EG{S)UEG{S)\^^_-^yS\(^^_-^-^nG-\S\ 
SQV{G) 

where the last equality is by considering the contribution of v' for each v e V{G): 
if V G S then v' contributes either — t or 1; if v ^ S then v' contributes either t or 
— 1. The last expression in the equation above equals 

(t-l)"° Yl tl-^°(^)^-^°(^)l = (t-l)"°-Z(G;t,l). 

SCV(G) 

□ 

3.2. The Ising polynomials of certain trees. We denote be Sn the star with n 
leaves. Let cent(6'„) be the central vertex of the star 5„. A construction based on 
stars will be used to interpolate the y indeterminate from Z{G\'y,5). First, notice 
the following: 

Proposition 20. For every n E N"*", 

Z(5„;{cent(5„)},0;t,y) = y • (yt + 1)" 

Z(5„;0,{cent(5„)};t,y) = (y + t)" 
Proof. By definition, 

Z(5„;{cent5„},0;t,y) = ^ ^\EsAS)uEsAsys\ 

S:{cent(Sn)}CSCV(S„) 

Z(5„;0,{cent54;t,y) = ^ ,\EsAS)uEsAS)\^\s\ 

SCy(S„)\{cent(S„)} 

Consider a leaf u of 5„. For Z{Sn; {cent(6'„)}, 0; t, y), a leaf v has two options: either 
V E S, in which case it contributes the weight of its incident edge, so its contribution 
is yt; 01 V ^ S, in which case it contributes 1. For 0, {cent(5„)}; t, y), v has 

two options: either v G S, in which case it does not contribute the weight of its 
edge, so its contribution is y; or u S, in which case its edge contributes t. □ 
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Figure 3. Examples of S-h and S-h[G). In (b), the black vertices 
u^v are the vertices of K2. They are also denoted cent('H,w) and 
cent('H,i') respectively. 

cent({2, 4}) cent({2, 4}, u) = u cent({2, 4}, v) = v 




(a) The graph S{2,4}- (b) The graph S'{2,4} (-f^2). 



Definition 21 (The graph S-u)- Let H be a set of positive integers. The graph 
S'H is obtained from the disjoint union of S'„ : n g "H and a new vertex cent(?^) by 
adding edges between cent('H) and the centers cent(S'„) of all the stars Sn ■ n £ H. 

See Figure [2Ia) for an example. 

Proposition 22. LefH he a set of positive integers. Then, 



Z(5«;{cent(H)},0;t,y) 

Z(5?,;0,{cent(H)};t,y) 

Proof. We have 

Z(5w;{cent(H)},0;t,y) 

Z(5„;0,{cent(H)};t,y) 



y-n (yt-(yt + i)" + (y + t)'0 



hen 

n 



y n 

hen 



n 

hen 



(y(yt + l)''+t.(y + t)'^) 

(t-Z(5,,;{cent(5,,)},0;t,y) 
+Z(5ft;0,{cent(5,,)};t,y)) 

(Z(^ft;{cent(^^)},0;t,y) 
+t-Z(^^;0,{cent(^,,)};t,y)) 



and by Proposition [501 the claim follows. 



□ 



Definition 23 (The graph S-h{G)). Let H he a set of positive integers and let G 
be a graph. For every vertex v of G, let S-H^y{G) be a new copy of S-h. We denote 
the center of each such copy of S-h by cent('H, v). Let S-h{G) be the graph obtained 
from the disjoint union of the graphs in the set 

{G} U {Sn,v ■■ V e V{G)} 
by identifying the pairs of vertices v and cent('H, v). 

In other words, Sy,{G) is the rooted product of G and (S'-^, cent('H)). See Figure 
for an example. 



Proposition 24. Let H be a set of positive integers. Let 



9p,n{'^,y) = 



n(y(yt + i)" + t.(y + t)'^)) 

.hen / 



\V{G) 



gy,n{t,y) = yJI 



hen 



yt ■ (yt + 1)^ + (y + t)'' 
y • (yt + 1)^ + t • (y + t)'' 
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Then 

Z(5„(G); t, y) = <7p.«(t, y) • Z{G- 1, 5y,w(t, y)) 
Proof. By definition 



Z(^«(G);t,y) = 

SCV(S-h(G)) 



^\Es^(G){S)UEs^iG){S)\y\S\ 



We would like to rewrite this sum as a sum over S C V{G). By the structure of 

Z(5«(G);t,y) = ^ ^\Ea[S)UEa{S)\ 
SCV(G) 



[]Z(^^,„;0,{cent(H,z;)};t,y) 



ves 



By Proposition!^ 

Z(5„(G);t,y) 



^ ^\EGiS)UEG{S)\ 

SCV{G) 



[(y n(yt-(yt + i)" + (y + t)")) 

\ \ hen ) 

( n(y(yt + i)'+t-(y+tf)) 



\S\ 



\V{G)\\S\ 



and the claim follows. 



□ 



The following propositions will be useful: 



Proposition 25. Let gy^'^{t.^y) be as in Proposition \24\ Let /iy,-H be the function 
given by 

^y,w(ei,e2,r) = J]^ ( 1 



hen 



ei + 62 • r" 



Let 7,(5 ^ {—1,0,1} such that 7 7^ —(5. There exist constants hi,ui,U2,w (which 
depend on 7 and S) such that for every two finite sets of positive even numbers Hi 
and 'H2 which satisfy 

• \Hi\ = andni,n2 C N+\{l,...,/ii}, 

we have 

(i) 9y,-Hi{l,S),gy^-Hi{l,S),hy^-Hi{ui,U2,w),hy^-H^{ui,U2,w) e M.\{0}, and 

(ii) 3y.-Hi(7,'5) = 5y,W2(7,'5) *if ^y,«i (^l, M2, w) = ^y.Wa (""i, ""2, w) 
Furthermore, ui and U2 are non-zero and w ^ { — 1, 0, 1}. 

Proof. It cannot hold that |(57 + 1| — \S + -fl. Furthermore we know that 7, (5 7^ 0. 
Hence, there is hi such that for every even h > /ii,the sequences S{S'y + l)^ + 7 • 
{S + 7)'' and ^7(^7 + 1)'' + {S + 7)'' are strictly ascending or descending, and in 
particular, are non-zero. Therefore we have ^y^-Hi (7, 5), 5y,-Hi (7, 5) G ]R.\ {0}. 
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We have for i — 1,2 



j-j <57'-(^7 + l)'+7-('5 + 7)'' 



11 5(^7 + 7. (,5 + 7)'' 



^ n 5(7^-1) •(57 + 1)'' 




(5((57 + l)'' + 7- ((5 + 7)'' 



hen ) 



Let ui — ^2"'_i ; W2 = and w = We have ui,U2,w G K \ {0} and 

w ^ {—1,1}. Hence, we can take hi to be large enough so that ui + U2 + w'' 
non-zero. Since ui +U2 + w^ is strictly ascending or descending for even h, we have 
^y.Wi (""ij W2, w), /iy,-H2 (wi, U2, w) G R \ {0} for large enough values of h. □ 

Proposition 26. Let j, 6 ^ { — 1,0,1} andj^~S. LetH be a set of positive even 
integers. Let (7p,?i(t,y)) be from Proposition \ 24\ Then there exists ft.2 such that if 
ncN+\{l,..'.,h2} then gp.nh, S)) ^ 0. 
Proof. Recall 

^ \ViG)\ 

5p,w(7,5)) ( 

We have that (5 + 7 is non-zero. If (^7 + 1 = then the claim holds even for h2 = 0. 
Otherwise, using that 1^7 + 1| ^ |i5 + 7I, at least one of {S^ + 1)'*, ((5 + 7)'* becomes 
strictly larger in absolute value than the other for large enough h. □ 

3.3. Proof of Theorem [Ij The following lemma is a variation of Lemma 4 in [B]. 

For any V., let cr(H) = J2heH ^■ 
Lemma 27. 

Let 7 ^ {-1,0, 1}, 5 7^ 0, 61,62 7^ and ri,r2 ^ {-1,0, 1} such that |ri| 7^ \r2\. 
For every positive integer q' there exist q = ^(q') sets of positive even integers 
'Ho, . . . , Jiq such that 

(i) a(H.) = 0(log3(z') for alii, 

(ii) a{Hi) = criHj) for all i ^ j, 

(iii) /t,-H,(ei,62,ri,r2) ^ /t,w^ (61, 62, ri, r2) for i ^ j, 

where /t,-H(6i, 62, ri, r2) is from Proposition [TS[ If additionally 5 ^ { — 1,1} and 
7 7^ —6, we have 

(iv) gy,n, {-1,5) ^ gy^H,{l,5) for i ^ j . 

(v) ffp,«,(7,'5)^0, 

where gy^u{eii ^2, ?'i) is from Proposition^^ 

The sets "Hi can be computed in polynomial time in q' . 

Proof. Let q — q' log^ q' . First we define sets Hq, . . . , H^. We will use these sets to 
define the desired sets Ho, ■ ■ ■ jHq- 

For i — 0,. . . ,q, let i[0], . . . , i[i] G {0, 1} be the binary expansion of i wher^ 
£ = [loggj. Let A be a positive even integer to be chosen later. Let r G {1,2} be 



■^Actually, we will also need that £ is larger than a constant depending on ei, but this is true 
for large enough values of q. 



COMPLEXITY OF ISING POLYNOMIALS 



13 



such that \rr\ = max{|ri|, 17-21}. Then |r3_^| — min{|ri|, |r2|}. Let mo be an even 
integer such that toq > hi from Proposition 1251 and toq > /i2 from Proposition 1261 
We choose as foUows: 

n[ = {mo + Afloggl • (2j + : < j < ^ ■ 

The sets "H- satisfy: 

a. they are distinct, 

b. they have equal cardinality 

c. they contain only positive even integers between mo 
and mo + /S.{\ogq + l)(21ogg + 1), and 

d. for i,j and any a £ and b £ Hp either a = 5 or 
\a-b\> Alogg. 

It is easy to see that cr{'H[) = n{\ogq), i = 0, . . . ,q. On the other hand, since all 
the numbers in each of the y.^ are bounded by 0(log^ q) and the size of each is 
O(logg), we get that cr{'H'^) — 0{\og'^ q) for each i. From this we get that at least 

q = Q, ^^T^T^^ = ^W) of the sets Hq, . . . ,1-1'^ have the same sum value cr('H-). 
Let {Hq^ . . . , 'Hq} be a subset of {Hq, . . . , "H^} such that all the sets in {Hq, . . . , Tiq} 
have the same sum value cr('Hi). We have (i), (ii) and (v) for "Hq, ■ • ■ ^Hq. 

We now turn to (iii) and (iv). The proofs of (iii) and (iv) are similar but not 
identical. 

Let < i 7^ j < g, 'Hi\j = "Hi \ Hj and 'Hj\i ~ Hj \ Hi- Notice n 'Hj\i — 0. 
Let a — a{T-Li\j) — and let d — \'Hi\j\ — \'Hj\i\. 

(iii) We write ft,-Hi for short instead of /t,-Hi (ei, 62, ri, r2) in this proof. When 
other parameters are used instead of ei, 62, ri, r2 we write them explicitly. 
Since /t,-H, — ft,u^\j ■ ft,-H,nHj, f\.,Hj = f'^Mj\i ' /t.W.nWj and Jt^-u^nUj 7^ 0, 
it is enough to show that /t.w,^^. - h,Hj\i 7^ 0- 
Since a{'Hi\j) — a{'Hj\i) we have 

/t,w.\, = /t,w,v iff /t,«,\^.(ei, 62, r^\rj^^) = /t,w,\. (ei, 62, rj^^) . 

Hence we can assume from now on that |r^| > 1 (otherwise we look at r^^ 
and instead). 

For every H, /t,-H can be rewritten as follows: 

hen XCH 

where 

.«(X)=(-^)'"'<W.3^i-\-). 

We think h E X (respectively h E T-L\ X) as corresponding to Crr!^ 
(respectively e^^rr^-r)- 
It suffices to show that 

(3) E E ^«A.w^o- 

It holds that SH^^^{n,\j) = sn^^,{n,\i) = Hence, 
s-Hi\j{T^i\j) and s-Hj-^i(Tij\i) cancel out in Inequality (jS]). Similarly, 
s-Hi\j{^) = sWj\i(0) = f^-T cancel out. Let mi be the minimal element in 
T-Li\j U T-Lj\i- Without loss of generality, assume mi e T-Li\j- We have 

^«,V,(HA. \ {^l}) = (^)'<^""^3"V • 
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SHi\j \ {'^i}) has the largest exponent of out of all the monomials 
in both of the sums in Inequality and any other exponent of is 
smaller by at least Alogg. We will show that Inequality ([3]) holds by 
showing the following: 

Each of the sums in Inequality (|H) has at most 2'°^*^+^ = 2q monomials 
corresponding to the subsets of 'Hi\j and 'Hj\i respectively. The abso- 
lute value of each of these monomials can be bounded from above by 
s ■ |r^|'^-™i-Aiog9|^g_^|mi+Aiog9^ ^j^^j.^ ^ ^j^g maximum of 

1. Hence, the right-hand side of Inequality ^ is at most 



and 



4q • s|r^ 



o"— mi — A log q 



r3- 



imi+Alogq 



63- 



Alogg 



in^\j - {mi})| . 



There exists a number A' which does not depend on q such that 
4^5 ^^j^^ < (A')'°s'^ and (iii) follows by setting A large enough so 



that A' • 



< 1. 



(iv) By Proposition [25l there exist ui, U2 7^ and w ^ { — 1, 0, 1} depending on 
7, S for which it is enough to show that 
hy^fi.{ui,U2,w) ^ /ly.-Hj (ui, U2, w) to get (iv). We write hy^y^. for short 
instead of /ly^-H; (wi, 1*2, w) in this proof. 

Since we have hy^-H^ = ^y,«.\, • ^y,«.n«j, ^y,«j = ^y.-Ha- ' ^y.«.n«j and 
hy.-HifTHj 7^ 0, it is enough to show that 

^y,Wi\j - /iy,-Hj\« ^ ^ ' 

i.e. 

1 



n 



1 



n 



1 



1 



Ul 4- U2 • w 



^0. 



or equivalently, 

(Ul+U2 - W ) ^ 

Consider a product of the form found in Inequality ([5])- 

J]^ (wi M2 • + 1) (Ul + U2 • W^) = 



XC-Hall-Ht, 



Let 



I^L„-(x) 
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It sufBccs to show that 



We have p(0) = p{'Hi\j U 'Hj\i) — 0, using that \'Hi\j \ — \'Hj\i\. Let mi 
be the minimal element in 'Hi\j U "Hj-y^ . Without loss of generality, assume 
mi e 'Hi\j- We have 



\p{Hi\jU'Hj\i- {■mi})\ = 
b({"''i})| = 



2(i— 1 2cr— m-i I 

9 (X/ 



|((wi + l)ui)'^-lu2ii''"H 



The largest exponent of w in Inequality ^ is -u;^"'^'"^ . For all other 
monomials in ([S]), the power of w is smaller by at least Alogg. Simil- 
arly, the smallest exponent of w in Inequality (O is w™^ . For all other 
monomials in ([6]), the power of w is larger by at least Alogg. 

Let Xq C 'Hi\j U 'Hj\i) be maximal with respect to |p(Xo)|. Since d < 
logq + 1, we can choose A large enough so that we have 
Xq = 'Hi\j U 'Hj\i — {mi\ if \w\ > 1 and Xq = {mi} if |ti;| < 1. 

We have the following: 



X<ZHi\jUHj\i:X^Xo 



< \p{Xo)\ 



implying that Inequality ^ holds. To see that Inequality ([7]) holds, note 
that 



E p(^) 



X^Xo 



< 2 



2 log q+2 



XlZH,\.jUHj\,:X=^Xo 



\p{x)\ 



Let 



k{d) = max + , l) • max (|wi| , l) • max ([uaj , l) 



Then 



E pi^) 



XCH^\jUHj\i: 
X^Xo 



< 



4 • 2l°S«+l • k{d) ■ |i„|2^-™i-Alogg^ > I 

4.2i°g9+i.fc(d) 



,|mi+A log q 



\w\ < 1 



So, there exists a constant c > depending on ui,U2,w such that (for 
large enough values of q), 



E 

XC■H,\jUHj^^,■. 

X^Xo 



P{X) 



< 



^log?|^|2£T-mi-Alogg^ |w| > 1 



\w\ < 1 



It remains to choose A large enough so that 

\w\ > 1 



A log g ^ {(ui+l)ui)'' 



\w\ < 1 



□ 
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We are now ready to prove Theorems [T] 

Theorem 28. Let (7,(5) € . If j ^ {-1,0, 1} and S^O, then 

(i) computing Z{G;"f,d) is -hard, and 

(ii) unless ^ETH fails, computing Z{G;"f,5) requires exponential time in 

"G 

log" no ■ 

Otherwise, Z{G; 7, 5) is polynomial-time computable. 

Proof. We set t = 7 and y = <5 with 7 ^ {— 1, 0, 1} and 6 ^i). By abuse of notation 
we refer to ci,C2,di,d2,\i, A2 from Lcmma[T3]as the values they obtain when t = 7. 
Since 7 7^ {—1,0,1}, it is easy to verify that the foUowing hold: 

a. ci+ di,C2 + d2 ^ 0, 

b. Ai,A2 ^0, 

c. Ai, A2 ^ ±(7^+7), and 

d. Ai ^ ±A2. 

Let Ci = '^'2j' and = ■y^^^ for i = 1,2. Let q' — Uq. Let Hq, • ■ • 1 "Wg be the 
sets guaranteed in Lemma B71 with respect to q' , 7, S, ei, 62, ri,r2. 

First we deal with that case that 7 7^ —5. We return to 7 = — (5 later. 

We want to compute the g + 1 values Z{G (g) Ht', 7, 1)- If 5 = 1 we simply do it 
using the oracle to Z at (7, 1). li S — —1 we use Lemma [TOl Otherwise we proceed 
as follows. 

By Proposition [Ml for each < i, k < q, 
(8) Z{Sn, {G®'Hk);i.5)= g^.n^ (7, S) ■ Z{G (^Ukn, gy,w. (7, S)) . 

It is guaranteed in Lemma [?7I that for i ^ j, gy -Hi(7,(5) ^ gy^-H^ (7, (5). 

We want to use Equation (O to interpolate, for each < fc < me, the univariate 
polynomials Z{G0T-Lk', 7, y)- We use the fact that the sizes of G^T-lk, and therefore 
the y-degrees of Z(G <Si Hk', "/,)/), are at most 0(nG log'^ nc). Since gp,Hi{l,S) is 
non-zero, we can interpolate in polynomial-time, for each < fc < mc the -f 1 
polynomials Z{G ig) 'Hfc; 7, y). 

So, we computed Z{G ^ 'Hk','J, 1) for < fc < q. Now we use these values to 
interpolate t and get the univariate polynomial Z{G (E) Hk, t, 1). By Lemma [T51 

Z{G; f,,n, (ei, 62, ri, r2), 1) = Z{G ® Hk; 7, 1) • ifp.u, (7))"' ■ 

Since 7 ^ {-1,0,1}, fp,-Hkh) ^ 0- By Lemma[271 /t,«fc (ei, 62, ri, r2) are distinct 
and polynomial time computable. Hence, the univariate polynomial Z{G;t, 1) can 
be interpolated. We get (i) by Proposition [51 Since Z{—;-f,S) is only queried on 
graphs S-H- (G ® Hk) of sizes at most 0{na log® no), (ii) holds by Proposition [TOl 
Consider the case 7 = —6. By Proposition [24l for every G we have 

Z(5{ij(G); 7, d) = {S-il- 6^)r^ ■ Z{G; 7, -5^) . 

and the desired hardness results follow by the corresponding for Z{G; 7, —S^) (using 
that 7 ^ -(-(5^), ^ {-1, 0, 1} and that {5 ■ (1 - J^))"^ ig non-zero). 

Now we consider the cases where 7 G {— 1,0,1} or 5 = 0. Two cases are easily 
computed, namely Z{G; l,S) = (1 (5)"'= and Z(G;7,0) = 1. 

The other two cases follows e.g. from Lemma 6.3 in [TJ]. In that lemma it 
is shown in particular that partition functions Za,d{G) with a matrix A of edge- 
weights and a diagonal matrix D of vertex weights can be computed in polynomial 
time if A has rank 1 or is bipartite with rank 2. For 7 = we have 
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Figure 4. The construction of the graph R^''^{P2) ior £ — 1 and 
g = 2, where P2 is the path with two vertices and one edge. 




so A is bipartite with rank 2. For 7 = — 1 we have 



1 -1 y \^ 1 

so A has rank 1. Note Lemma 6.3 extends to negative values of S. We refer the 
reader to [T3] for details. □ 

4. Simple Bipartite Planar Graphs 

In this section we show that the evaluations of Z{G; x, y, z) are generally #P-hard 
to compute, even when restricted to simple graphs which are both bipartite and 
planar. To do so, we use that for 3-regular graphs, Z(G';x, y,z) is essentially equi- 
valent to Z{G; t, y). We use a two-dimensional graph transformation R^''^{G) which 
in applied to simple 3-regular bipartite planar graphs and emits simple bipartite 
planar graphs in order to interpolate Z(G;t,y). 

4.1. Definitions. The following is a variation of k-thickening for simple graphs: 

Definition 29 (fc-Simple Thickening). Given £ S N+ and a graph iJ, we define 
a graph STh^{H) as follows. For every edge e = {u,w) in E{H), we add M new 
vertices Ue,i, . ■ . , We,« to H. For each Ve,i, we add two new edges {u,Ve,i) and 
(WyVg^i). Finally, we remove the edge e from the graph. Let Ne{e)~^ denote the 
subgraph of STh^(H) induced by the set of vertices {we,!, • ■ • , We,«, u, v}. 

The graph transformation used in the hardness proof is the following: 

Definition 30 (i?^^'?(G)). Let G be a graph. For each w e V{G), let = (V^, E^) 
be a new copy of the star with 2q leaves. Denote by Cw the center of the star G^. Let 
= {V^^e'jI'') be the graph obtained from the disjoint union of STh^{G) 
and STh\Gl,) for all w G V{G) by identifying w and c„ for all w € V{G). 

Remarks 31. 

(i) The construction of R^'''{G) can also be described as follows. Given G, we 
attach 2q new vertices to each vertex v of V{G) to obtain a new simple 
graph G' . Then, R^^'^iG) = STh^{G'). 

(ii) For every simple planar graph G and i,q E N^, R^''^{G) is a simple 
bipartite planar bipartite graph with nj^ vertices and rriR edges, where 
nR — nc{l + 2g(l + 4£)) + A£ma and and mji — 8imc + IGiqno- 

Figure m shows the graph R^'^{P2). 

In the following it is convenient to consider a multivariate version of Z(G] x, y, z) 
denoted Z{G; x, y.z). This approach was introduced for the Tutte polynomial by A. 
Sokal, see [13]. Z{G; x, y.z) has indeterminates which correspond to every v € V{G) 
and every e € E{G). 
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Definition 32. Let x {xe : e e E{G)), y = (y„ : u e F(G)) and 

z = (ze : e e E{G)) be tuples of distinct indeterminatcs. Let 

Z(G;x,y,z)= E ( 11 flly") ( 11 ' 

5Cy(G) \eeEG[S) J \ueS I \eeEG{S) ) 

We may write x^ and ~Lw,v instead of Xg and Zg for an edge e — {w, v). Clearly, 
by setting Xg = x and — z for every e G E{G), and = y for every u G V{G) we 
get Z(G;x,y,z) = Z(G;x,y,z). 

We furthermore define a variation of Z(G;x,y, z) obtained by restricting the 
range of the summation variable as follows: 

Definition 33. Given a graph H and B.C C V{H) with B and G disjoint, let 
(9) ZiH,B,C;x,y,z)^ 

E f n n y-l f n ^] 

A:BCACViH) ,AnC=(/) \eeEG(A) j \ueA\B j \e.eEa{A) ) 

where the summation is over all A C V{11\ such that A contains B and is disjoint 
from G. 

We have Z{H, 0, 0; x, y, z) = Z{H; x, y, z). 

4.2. Lemmas, statement of Theorem [2] and its proof. For every edge e G 
E{G) between u and v, let 

uJi{e,S) = Z(7VKe)+,S'n{u,i;},{M,w}\S';x,y,z), 

and for every vertex w G V, let 

uj2{w, S) ^ Z{STh\Gl), S n M, M \ S; X, y, z) . 

Let 

^il-S*) = H ^^1(6, 5") and uj2iS) ^ u!2{w,S). 

eeE(G) w£V(G) 

Let Wi^trivl'S') for i = 1, 2 be the polynomials in x, y and z obtained from u}i{S) by 
setting Xe = X and z^ = z for every e G E^'*^ and y„ = y for every u e V^'' . 

Lemma 34. 

Z(i?^^''(G); x.y.z) = ^ c.i,triv(5) • W2,triv(^) • yl^l . 
Scy(G) 

Proof. Each edge of R^''^{G) is either contained in some A'£(e)+ for e G E{G), or 
in some STh^{Gl) for w e V^(G). Hence, by the definitions of Z(i?^'«(G); x, y, z), 
wi(S') and ^2(5'), 

Z(i?^'«(G);x,y,z)== ^ uj,{S) ■ Lo^iS) ■ \{y^ 
Scy(G) toes 
holds and the lemma follows. □ 

Lemma 35. Let e = {u,w) be an edge of G. Then 

r(y + z2)« \{u,v}nS\=0 
wi,triv(e, S) = Uxy + z)« \{u, v} n S\ ^ 1 
[(yx2 + l)« \{u,v}nS\^2 

Proof. The value of wi(e, S*) depends only on whether u,w G S. Consider A C 
V{Ni{e)~^) which satisfies the summation conditions in Equation ^ for Z{Ni{e)^,Sr\ 
{u,w},{u,w}\S;x,y,z). 
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(i) li w £ S and u ^ S: Exactly one edge e' incident to Ve,i crosses the cut 
[A, A]^^(e)+- The other edge e" incident to We,i belongs to E{A) or E{A), 
depending on whether Ve,i S A. We get: 

Wi(e, S) = JJ(x^,,i,u;y!;,,i + Zt>e.,,«) ■ 
i=l 

(ii) If w ^ S* and u G S: this case is similar to the previous case, and we get 

i=l 

(iii) If w, It € S': For each We_i, either v^j G A, in which case both edges (we.i, w) 
and (we,i,u) are in E{A), or We,i ^ 5, and both edges {ve,i,'w) and {ve.i,u) 
cross the cut. We get: 

4£ 

u}i{e,S) = Y]_{Vv,.i><v^,,,uy-v^,,,w + 1) • 

i=l 

(iv) If u ^ 5': For each Ue.i, either v^j £ S and then both edges incident to 
Ve,i cross the cut, or v^.i ^ S and none of the two edges cross the cut. We 
get: 

i=l 

The lemma follows by setting Xe = x and Zg = z for every edge e and Yu — y for 
every vertex u. □ 

Lemma 36. Let 

ff£,g(x,y,z) = y-(yx +1) +(yx + z) 
^f,g(x,y,z) = (y + z^j^Vy- (yx + z)"*^ . 
Let w be a vertex of G. Then 



W2,triv(w,S') = 



(.9^9(x,y,z))^'' w£S 
{h,^,{x,y,z)f' wis 



Proof. Consider A which satisfies the summation conditions in Equation ^ for 
Z{STh\Gl), S n {w}, {w} \ S- X, y, z). 

(i) If € S* (or, equivalently, c^, € A): Let u £V^ \ {c^,} and e — {w, Cm}. If 
u € A, then the vertices u and We,i, • . • , We,4£ contribute 



41 



Yu J_ J_ (y^'ci Xtfg^^ ^W^Ve,i ,U ^" I) ■ 

i=l 

Otherwise, if u ^ A, then the vertices u and Wg^i, . . . , contribute 

4£ 
i=l 

Hence, uj2{w, S) equals in this case 

C4£ 4e \ 

+ 
2-1 2-1 / 
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(ii) M w ^ S (or, equivalently, ^ A): Let u G V^\ {cw} and e = {u, Cw}- If 
u G A, then the vertices u and We,i, . . . , t;e,4£ contribute 

4i 

1 = 1 

Otherwise, if u ^ A, then the vertices u and Ve.ii ■ • ■ , We,« contribute 

i=l 

Hence, lo2{w, S) equals in this case 

C4e 4£ \ 

The lemma follows by setting Xg = x and Zg = z for every edge e and Yu = Y for 
every vertex u. □ 

Lemma 37. If G is d— regular, then 

fpM^, y, z, e, q) ■ Z[G- /t,fl(x, y, z, ^), /y,fl(x, y, z, g)) = Z(i?^^«(G); x, y, z) 

w/iere 

/p,«(x,y,z,£,q) - (/i,,,(x,y,z))^«"«(y + z2)2"«o^ 

f ^ «^ ( (yx + z) 

/t,i?,(x,Y,z,£) = 



X 2f 



/y,-R(x,y,z,£,q) = y 



(yx2 + l)(y + z2) 

yx + 1 A / g<;,g(x,y,z) ^ 
y + z2 y V/i^,,(x,y,z) 



Proof. We want to rewrite Z{R^'''{G); x, y, z) as a sum over subsets 5* of vertices of 
G. Using Lemma [Ml in order to compute Z(i?^'^(G); x, y, z) we first need to find 
wi,triv('5') and W2,triv('5')- Using Lemma [36l CL;2,triv(S') is given by 

^2,triv(5) = (g,.,(x,Y,z))2'l^l . (/i,,,(x,Y,z))2^"^-2«l^l . 

In order to compute a;i^triv(5'), consider 5* C V{G). Since G is d-regular, the number 
of edges contained in 5 is ^{d- \S\ — \ [S,S]g\), and the number of edges contained 
in S is ^{dnc — d - \S\ — \ [S, S]g\)- Hence, by Lemma [35l wi^triv(5') is given by 

^i.tMS) = (xy + z)«l[^'^l-l (yx^ + 1^ ""'"'^'""' (y + ,2)«.±i.i^±^HiMM _ 

Using Lemma [Ml 

Z{R'''^{G)- X, y, Z) = ^l.triv(^) • W2,triv(^) ' yl'^l 

SCV{G) 

which is equal to (y + z^)^^ ~^ times 

/ rvx + 7^2 X 2£|[S.S]gI / /„^2 , TX2M\ 1-51 

Plugging the expression for a;2,triv(5') in Equation ((TU]), we get that Z{R^''^{G); x, y, z) 
equals /p,_r(x, y, z, g) times 



E 

SCV{G) 



I , .2 X 2£|[S,S]gI / / 2 , 1 \ 2^d . I s X 2g\ 

(yx + z)^ \ ^ '/ /yx^ + l\ /g£,,(x,y,z) ^ 



,(yx2 + i)(y + z2); Vy + zV v/i«,g(x.y.z), 

and the lemma follows. □ 
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Lemma 38. Let e E Q\{— 1,0,1} and let a,b,c > and b ^ c. Then the 
ci e N for which the sequence 

(11) m ^ 



+ e • a* 

is strictly monotone increasing or decreasing for £ > Ci. 
Proof. h(£) can be rewritten as 

e-b^ + 1 

by dividing both the numerator and the denominator of the right-hand side of 
Equation (fTT|) by and setting b = ^ and c = ^. We have b ^ c and b,c> 0. 

Let h{x) = "^gx ■ The derivative of h{x) is given by 

e In 6 • 6^ (P + e) — In c • P (^e ■b'^ + 1^ 



h'{x) 



hi 6 • &^ — In c • P + e(ln 6 — In cjb'^c^ 



The denominator of h'{x) is non-zero for large enough x. Therefore, there exists xq 
such that h'{x) is continuous on [xq, oo), so it is enough to show that h'{x) ^ for 
all large enough x to get the desired result. 

If 6 = 1 then + e)^ h' (x) = -(1 e) InS • P, and if S = 1 then (2^ + e)2/i'(a;) = 
(e^ + e) ln6 • 6^. In both cases h'{x) is non-zero, using that 6^5 and b,c > 0. 

Otherwise, 6, c and be are distinct. Let Ai = {b^,c^,b^c^}. Let A2 be the 
subset of Ai which contains the functions of Ai which have non-zero coefficients in 
Equation Note belongs of A2. There is a function in A2 which dominates 
the other functions of ^2- This implies that h'{x) is non-zero for large enough values 
of x. □ 

Theorem [5] is now given precisely and proved: 

Theorem 39. For all (7, S, e) G such that 

(i) J ^{-1,0,1}, 

(ii) J + ^{-1,0,1}, 

(iii) (5 + eV±(^7' + l), 

(iv) -572 -f 1 =^ 0, 

(v) j6 + e ^ 0, and 

(vi) (75 + e)V ('57' + l)"!-^ + 

Z(— ;7, (5, e) is ^V-hard on simple bipartite planar graphs. 

Proof. We will show that, on 3— regular bipartite planar graphs G, the polynomial 
Z(G\ t, y) is polynomial-time computable using oracle calls to Z(— ; 7, 5, e). The or- 
acle is only queried with input of simple bipartite planar graphs. Using Proposition 
[HI computing Z(G';t,y) is #P-hard on 3-regular bipartite planar graphs. 

Using ^ and ^ it can be verified that there exists cq € N+ such that for all 
£>Cq and q € N+, fp,R{"l, 5, e, 2f , q) ^ 0. We can use Lemma [571 to manufacture, in 
polynomial-time, evaluations of Z{G; t, y) that will be used to interpolate Z{G; t, y). 

Let £> Co and let 

We have that fy^ji{j,5,£,l,q) = (5(_Ey4)'''^ (Ey^2.i)^''- Using ([Iv]) we have i?y,i 7^ 0. 
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Look at Ey_2,g as a function oi I. Using (P, ((n]), (pli)) and (jivl) and Lemma [551 with 
a = + 6 = ((57^ + 1)'*, c = + and e = i5, there exists ci such that £'y,2,^ 
is strictly monotone increasing or decreasing. Hence, there exists C2 > ci such that, 
for every ^ > C2, -Ey,2,£ ^ {—1: Oj !}• Moreover, C2 = 02(7, (5, e) is a function of 7, (5 
and e. 

We get that for gi 7^ 92 € ["g + 1] and t > ca, (i^y,2,f)^'' 7^ (£^y,2/)^'''. Since 
S {Ey^if''^^ is not equal to and does not depend on g, we get that for qi 7^ (72 S 
[no + 1], fy,Bi7,S,e,e,qi) ^ /y,fl(7, <5, e, ^, 52)- 

For every i G [ma + C2 + 1] \ [02], wc can interpolate in polynomial-time the 
univariate polynomial Z(G; /t,fl(7, (5, e, ^), y). Then, we can use the polynomial 
Z{G; ft,R{^,S,e,e),y) to compute Z{G; ft,R{'j,S,eJ),j) for every i e [mc + C2 + 
1] \ [C2] and every j £ [no + !]■ Let 

^ / (7^+s)^ y 

* V('57' + l)('^ + e')/ 
and it holds that ft,R{j,5,e,i,q) = {E^Y . Clearly, Et ^ —1 and, by (jvj) and (fvil) . 

^ {0, 1}. Hence, for every ^ £2 e N+ we have /t,fl(7, ^, e, ^1) 7^ /t,_R,(7, e, ■^2)- 
Therefore, we can compute the value of the bivariate polynomial Z(G; t, y) on a grid 
of points of size (tog + 1) x (ng + 1) in polynomial-time using the oracle, and use 
them to interpolate Z(G;t,y). 

□ 

5. Computation on Graphs of Bounded Clique-width 

In this section we prove Theorem[3] Let G be a graph and let cw{G) be its clique- 
width. As discussed in Section^ a ^-expression t{G) for G with k < 23^™('^)+2 - 1 
can be computed in FPT-time. Let c = (c„ : w € ^{G)) be the labels from 
[k] associated with the vertices of G by t(G). We will show how to compute a 
multivariate polynomial Ziabcicd(G, c; x, y, z) with indeterminate set 

{x{-*,j},yi,z{,,j} I i,j e [k]} 

to be defined below. Note it is not the same multivariate polynomial as in Section 
21 For simplicity of notation we write e.g. Xi,j or Xj^i for x^ij} ■ The multivariate 
polynomial Ziabcicd(G, c; x, y, z) is defined as 

(13) E fny^^)f n xe^.^U n ^^^^X 

The left-most product in Equation (|13p is over all vertices v in S. The two other 
products are over all edges in Ec{S) and Ea{S) respectively. It is not hard to 
see that Z(G;x, y,z) is obtained from Ziabeied(G, c; x, y, z) by substituting all the 
indeterminates Xij, y^ and z^j by three indeterminates, x, y and z, respectively. 

Given tuples of natural numbers a = {oi : i G [k]), 6 = (4.j : hi ^ [k]) and 
c = {ci,j : «, j G [fc]), we denote by t-^j^iG) the coefficient of the monomial 

n yf n <Mj 

in Ziabciod(G; X, y, z). We call a triple (a, 6,t) valid if ai + . . . + o/j < uq and, for 
all i,i £ [k], 6ij,Cij < niQ. If {a,6,c) is not valid, then t--^-{G) = 0. Therefore, 
to determine the polynomial ^iabcicd(G; x, y, z) we need only to find i j(G) for all 
valid triples (a, 6, c). 

The t--i^-{G) form an (A; 2fc^)-dimcnsional array with (max {rtc, tog})'^^^'^ 
integer entries. Each entry in this table can be bounded from above by 2""^ and 
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thus can be written in polynomial space, so the size of the table is of the form 
^g,pi(cw(G))^ where pi is a function of c'w{G) which does not depend on no- 

We compute 2'iabciod(G, c; x, y, z) of G by dynamic programming on the structure 
of the fc-expression of G. 

Algorithm 40. 

(i) // (G, i) is a singleton of any color i, 

■^labeled 

(G,c;x,y,z) = 1 + yj. 

(ii) If{G,c) is the disjoint union of (H, chi) and {H2,ch2), then 

^labeled (G, c; X, y, z) — Ziaboicd(^^i, c//i ; X, y, z) • Ziaboicd(-ff2,c_ff2;x, y, z) . 

(iii) // (G, c) = rjp_r{H, ch): let cCr and dp he the number of vertices of colors r 
and p in H , respectively. 

(iii. a) For every valid (a, 6,1), if 

j Op- Or p=^r { {tip - Op) ■ {dr - Or) p ^ r 

(14) ^--|(^) '^"^ ''^'^\{'^-,^) V = r 



set 



where the summation is over all valid tuples 6' — {6[ ^ : i,j & [k]) 
and ~c! = [c[ ^ : i,j € [k]) such that 6[ ^ = 6i,j and c- ^ = c^j if 
{hj}¥'{P,r}. 

(iii.b) For every valid (a, 6, t), if Equation does not hold, set t-^ i j(G) = 
0. 

If{G,c)=pp^r{H,CH): 

(iv.a) For every valid (a, 6, c) if Op = 0, set 



a' ,B' .c' 

where the summation is over all valid tuples a' = {al : i € [k]), 
^' = ■ hi ^ [k]) and c' = {clj : i,j £ [k]) such that 

• Or = Up + fl'.. 

• CH = a- for all i ^ {p, r}, 

• for all j e [k] \ {p}, 




— ) 3,P ' 3,r 



C. 



'r,r ~^ ^p,r ~^ ^p,p J ^ 



and 

• for all i,j £ [k] \ {p,r}, = and — c'^ ^ 
(iv.b) For every For every valid (a, 6, c) if Op ^ 0. set tj ^ j(G) = 0. 

Correctness: 

(i) By direct computation. 

(ii) Proved in [1 for Z(G;t,y). The trivariate case is similar. 

(iii) G = rjp r{H): Let 5 be a subset of vertices of V{G) — V{H) with 
Op and Or vertices of colors p and r respectively. After adding all 
possible edges between vertices of color p and of color r in S", the 
number of edges between such vertices in Eg{S) is Op ■ Or ii r p 
and ('^) if p = r. Similarly, the number of edges between vertices 
colored p and r in EaiS) is {dp — Op) ■ {dr — Or) ii r ^ p and (''^ 2 

if p = r. 



24 



TOMER, KOTEK 



(iv) G = Pp^r{H): Let be a subset of vertices of V{G) = V{H). 
After recoloring every vertex of color p in to color r, we have 
Op = 0. Every edge between a vertex colored p to any other 
vertex lies after the recoloring between a vertex colored r and 
another vertex. There is one special case, which is the edges 
that lie between vertices colored r after the recoloring. Before 
the recoloring these edges were incident to vertices colored any 
combination of p and r. 
Running Time: The size of the (^2^<:'"(G)+2 _ i)_expression is bounded by 
n'^ ■ fi{k) for some constant c, which does not depend on cw(G), and 
for some function /i of cw{G). Now we look at the possible operations 
performed by Algorithm HUl 

(i) The time does not depend on n since G is of size 0(1). 

(ii) The time can be bounded by the size of the table ^ ^ to the 
power of 3, i.e. n^Pi(^-^iG)) , 

(iii) For fip^r, the algorithm loops over all the values in the table tj ^ j, 
and for each entry possibly compute a sum over at most mc ele- 
ments. Then, the algorithm loops over all the values again and 
performs 0(1) operations. 

(iv) For pp^r, the algorithm loops over all the values in the table ^ j, 
and for each entry possibly compute a sum over elements of the 
table t--^-. Then, the algorithm loops over all the values again 
and performs 0(1) operations. 

Hence, Algorithm HUl runs in time O (^Hq™'"'^^''^ for some function /. 

6. Conclusion and Open Problems 

Applying the reductions used in the proof of Theorem [T] to planar graphs gives 
again planar graphs. Combining Theorem [T] and its proof with Lemma [571 a hard- 
ness result for the trivariate Ising polynomial on planar graphs analogous to The- 
orem [5] follows. However, both Theorem [5] and the analog for planar graphs are 
not dichotomy theorems since each of them leaves an exceptional set of low dimen- 
sion unresolved. Theorem [2] serves mainly to suggest the existence of a dichotomy 
theorem for Z(G;x, y,z) on bipartite planar graphs. 

Another open problem which arises from the paper is whether Z(G;x, y,z) re- 
quires exponential time to compute in general under ^ETH. One approach to the 
latter problem would be to prove that, say, the permanent or the number of max- 
imum cuts require exponential time under ^ETH even when restricted to regular 
graphs. 
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